function y = liarwhd(x)
% LIARWHD function
% x0=[4,4,...,4];
n=length(x);
y=0;tp1=0; tp2=0;
for i=1:n
    tp1=tp1+4*(-x(1)+(x(i))^2)^2;
end
for i=1:n
    tp2=tp2+(x(i)-1)^2;
end
y=y+tp1+tp2;


